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We present an analytical model for describing highly relativistic acceleration 
of magnetically driven jets, within the framework of ideal MHD for cold, sta- 
tionary and axisymmetric outflows. Our novel procedure is to treat the wind 
equation as an algebraic relation between the relativistic Alfven Mach-number 
and the poloidal electric to toroidal magnetic field amplitudes ratio £. This allows 
us to obtain easily the wind solutions for trans-fast-magnetosonic flows, together 
with the required range of £. Then, to determine the spatial variation of £, we 
solve approximately the Grad-Shafranov equation applied to a jet flow ejected 
with a very large total specific energy E and confined within a very small opening 
angle. Our trans-fast-magnetosonic model provides a closed-form expression for 
the transition from a magnetically dominated flow to a kinetic-energy dominated 
one, which occurs in the sub-asymptotic region far beyond the light cylinder of 
the radius R^. Importantly, we find that the equipartition between magnetic 
and kinetic energies is realized at a cylindrical radius of order of R^E/c 2 , and 
confirm that the further conversion of magnetic to kinetic energy proceeds loga- 
rithmically with distance in the asymptotic region. Finally, we discuss briefly the 
astrophysical implications of our model for jets originating from active galactic 
nuclei. 
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1. INTRODUCTION 

1.1. Collimated Relativistic Jets 

Highly relativistic jet motion is one of the most interesting phenomena observed in ac- 
tive galactic nuclei (AGN), microquasars and possibly 7-ray bursts (GRB), and many works 
have been devoted to numerical and analytical investigations of the mechanisms for produc- 
ing, collimating and accelerating matter to relativistic speeds (Sauty et al 2001; Ghisellini 
& Celotti 2002). At present the most promising approach to the jet phenomena seems to 
consider magnetically driven outflows within the framework of relativistic magnetohydrody- 
namics (MHD). 

Based on the MHD scenario, the jet ejection is expected to be realized under the mag- 
netic energy-dominated state (i.e., the Poynting jet) in the vicinity of a central source. This 
initial Poynting flux will be originated by electromagnetic extraction of rotational energy 
from a spinning black hole or/and an accretion disk, as was first discovered by Blandford 
& Znajek (1977) and numerically studied by the time-dependent MHD simulations in Kerr 
geometry (see Komissarov 2001 for the magnetically dominated regime; Koide et al. 2002 
for the full MHD regime). Then, a significant fraction of the huge energy in the outflow will 
be converted from the Poynting flux into the fluid kinetic energy of bulk motion. Such an 
efficient energy conversion may be able to occur unsteadily in the inner magnetosphere close 
to the black hole as a result of the MHD interaction with infalling matter. Koide et al. (2000) 
discussed this problem and found a magnetically driven jet inside a gas pressure-driven jet in 
the counter-rotating black hole case against the disk rotation. However, the poloidal velocity 
of the jet is only sub-relativistic. 

Otherwise, the steady magneto-centrifugal acceleration in the propagation over a large 
enough distance should become important. It is well-known that the ideal stationary axisym- 
metric MHD equations reduce to a set of two equations describing the local force-balance 
along the field and across the field, and called the poloidal wind equation and the Grad- 
Shafranov (trans-field) equation (for reference, see, e.g., Okamoto 1992; Beskin 1997). Self- 
similar solutions of magnetosphere, where the magnetic field lines are anchored to a thin 
accretion disk, were discussed for jet collimation and acceleration (Blandford & Payne 1986; 
Li, Chiueh & Begelman 1992; Contopoulos 1994). Unfortunately, the outflows given by 
the self-similar solutions (Contopoulos 1994) reach some maximum radius and recollimate 
afterwards. The self-similar scaling will not be valid in the asymptotic region. Using the 
stationary axisymmetric model, the asymptotic flow structure has been found to vary log- 
arithmically with radius, and the kinetic energy- dominated solutions describing collimating 
jet magnetospheres have been presented (Li, Chiueh & Begelman 1992; Eichler 1993; Begel- 
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man & Li 1994; Tomimatsu 1994), while the decollimation of magnetic field lines has been 
claimed to become significant in the equatorial (current-closure) region (Beskin & Okamoto 
2000). Note that the asymptotic analysis fails to check whether the obtained structure is 
due to trans-magnetosonic flows satisfying the critical conditions. To assure the validity of 
the asymptotic structure, we should consider the connection of the plasma source to the 
asymptotic region (and to the black hole). However, the full MHD description of the global 
magnetospheric structure in general relativity is a task still too complex even for the sta- 
tionary axisymmetric system. For example, the available general relativistic models to give 
both closed and open field lines connecting a disk with a black hole and a distant region 
are limited to magnetic-energy dominated cases, such as a disk-current field (Tomimatsu & 
Takahashi 2001) and a force-free field (Fendt 1997). 

The jet-type outflows given by the one-dimensional or two-dimensional numerical non- 
self-similar solutions, on the other hand, were also constructed (Appl & Camenzind 1993a,b; 
Li 1993; Fendt & Camenzind 1996; Fendt & Memola 2001). By solving the Grad-Shafranov 
and poloidal equations self-consistently, disk-jet connected flow solutions cylindrically con- 
fined by some external pressure were obtained. To compare the MHD predictions with jet 
observations in more detail ( e.g., the terminal velocity and width of the jets, the size of the 
active disk region in a central engine, which would correspond to jet's foot point, and the 
gravitational radius ), the efficiency of relativistic acceleration should be discussed in more 
general asymptotic structures of outflows (depending on the boundary conditions; e.g., the 
shape of the innermost flux surface, or an external pressure distribution at the outermost 
flux surface, etc), and the propagation distance necessary for a rough equipartition between 
magnetic and kinetic energies should be clarified. In this paper we would like to focus on 
such problems of steady acceleration at large distances from the central source, by developing 
a new analytical method without any self-similar assumption. Most of the previous works 
discussed the field structure under the force-free assumption, in which the contributions by 
the plasma inertia are ignored. However, our new treatment allows us to study more easily 
trans-fast magnetosonic outflow solutions and to include the plasma inertia effects on the 
magnetic field structure and the plasma acceleration process. Hence, by neglecting the effect 
of gravity and gas pressure, we can have some more self-consistent MHD solutions, showing 
the evolution of the bulk Lorentz factor 7 of outgoing flows to 7 ~ 10 in accordance with 
AGN jets. 
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1.2. New Analytic Procedure for Trans-Magnetosonic Jets 

In this paper, for the ideal stationary axisymmetric MHD outflows, we follow the rela- 
tivistic formulation given by Camenzind (1986, 1987). To avoid a very troublesome task to 
solve the Grad-Shafranov equation, one may discuss the fluid motion by assuming a mag- 
netic field configuration and solving only the poloidal wind equation (Takahashi & Shibata 
1998; Fendt & Greiner 2001). Though this procedure does not always assure the consistency 
with the Grad-Shafranov equation, the analysis of the poloidal wind equation is useful for 
studying the critical conditions for MHD flows. For a cold plasma considered here, the key 
problem is to derive the trans-fast-magnetosonic solutions from the poloidal wind equation. 
If the poloidal wind equation is algebraically solved under a fixed field configuration in a con- 
ventional manner, however, it allows the existence of many unphysical solutions describing a 
flow which does not smoothly pass through the fast-magnetosonic critical point, unless the 
slightly complicated critical conditions are required for the four integrals of motion which are 
constant along a filed line. Usually this prevents us from understanding clearly the efficiency 
of acceleration along a field line in simple analytic computations. 

One of our purposes in this paper is to propose a new analytic procedure useful for 
discussing general properties of the trans-MHD flows within the poloidal wind equation and 
without assuming a field configuration. We treat the poloidal wind equation as an algebraic 
relation between the relativistic Alfven Mach-number M and the poloidal electric to toroidal 
magnetic field amplitudes ratio £ (instead of the poloidal field amplitude). Though it is not 
difficult to contain the effect of gravity, our analysis given here is limited to special relativistic 
flows for showing explicitly the conversion of magnetic to kinetic energy at large distances. 
Then, it is easily found that only if £ is supposed to be a function smooth along a field line 
at the critical points, the solution (depending on £) for the poloidal wind equation becomes 
trans-critical. We would like to emphasize that no special critical conditions are needed for 
the four integrals of motion. 

The spatial variation of the ratio £ newly introduced into the poloidal wind equation 
must be determined by the Grad-Shafranov equation. In the case of AGN (or possibly GRB) 
jets a high degree of collimation as well as a ultrarelativistic bulk motion is observationally 
indicated. Hence, we consider jet flows with a very small opening angle and a very large 
asymptotic Lorentz factor to solve approximately the Grad-Shafranov equation, under the 
condition of a cold plasma and no gravity. Hereafter we use the cylindrical coordinates R 
and Z with c = 1 unit and denote the flux function by ty(R,Z). For efficient acceleration 
of bulk motion a huge magnetic energy should be stored in the outflow near the central 
source. Then, we assume the total specific energy E(^f) to be much larger than unity. (In 
this paper, we use the term "specific" as a quantity per the rest-mass energy of one plasma 
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particle.) The specific energy Efif) is conserved along a field line ^ = constant. The 
angular velocity of a field line is also an integral of motion, and the light cylinder 

radius is given by = 1/Vtp. By virtue of the very large value of E, we can define the 
intermediate region <C R < R\^E between the inner region R < R^ and the asymptotic 
region R ^> R^E. We derive a class of general solutions corresponding to an arbitrary (i.e., 
cylindrical, paraboloidal or radial) boundary configuration for the outermost flux surface 
of jet flows with the opening angle smaller than 1/E, by solving the approximated Grad- 
Shafranov equation, which remains valid in the spatial range from the intermediate region 
to the asymptotic region. The total specific energy E may be decomposed into the magnetic 
part E m and the plasma kinetic part E k . The importance of the energy conversion from 
the magnetic part to the kinetic part in a relativistic wind was originally pointed out by 
Michel (1969), and he showed that for a radial magnetic field the energy conversion from 
Poynting flux to particle energy flux (or the plasma acceleration) is extremely inefficient. 
For a flux tube diverging even slightly faster than radially, the significant energy conversion 
occur beyond the fast magnetosonic point; at the fast magnetosonic point the total energy 
flux is still dominated by the Poynting flux (e.g., Begelman & Li (1994); Takahashi & Shibata 
(1998)). The main purpose of this paper is to give the evolution of the ratio E k /E m in a 
simple closed-form expression, and our new result is to show the conversion of E m to E k 
in the intermediate region. The fast-magnetosonic critical point is found to exist in the 
intermediate region, where the flow is still magnetic-dominated (E m ^> E^) in accordance 
with the previous works. Importantly, we can conclude that the rough equipartition E m ~ Ej, 
is realized at a radius R ~ RlE. The further energy conversion goes on logarithmically with 
the increase of R toward the kinetic-energy dominated asymptotic state, as was previously 
pointed out (Eichler 1993; Begelman & Li 1994; Tomimatsu 1994). The complete energy 
conversion becomes possible, if the outermost flux surface of jet flows extends to an infinite 
radius. 

This paper is organized as follows: In §2, we analyze the poloidal wind equation by 
introducing the poloidal electric to toroidal magnetic field amplitudes ratio £. The general 
properties of the trans-fast-magnetosonic solutions are studied, in particular, for outflows 
which are initially ejected under a magnetic-energy dominated state and a very large total 
specific energy E. In §3, we give the approximated form of the Grad-Shafranov equation for 
jet flows propagating beyond the light cylinder surface with a very small opening angle. In 
§4, we derive a class of general solutions for the Grad-Shafranov equation to determine the 
Alfven Mach number M and the flux function \l/ in the intermediate and asymptotic regions. 
Finally we summarize our results of relativistic acceleration and discuss the implications for 
AGN jets in §5. 



2. TRANS-CRITICAL MHD OUTFLOWS 



2.1. Non-singular MHD flow equation 

First let us give a brief review (see, e.g., Camenzind 1986) of ideal stationary axisymmet- 
ric MHD outflows for a cold plasma in Minkowski spacetime background with the cylindrical 
line element written by c = 1 unit as follows, 

ds 2 = dt 2 - dR 2 - dZ 2 - R 2 d(f) 2 . (1) 

As was mentioned in §1, we have the four integrals of motion which are constant along a field 
line given by ^(R,Z) = constant. The total specific energy E(^f) and angular momentum 
L(\E') are decomposed into the magnetic and kinetic parts as follows, 

E = E m + E k , E m = ^ - , E k — 7 , (2) 

RB 

L — L m + L k , L m = — - — — , L k = 7/? , (3) 

Auk 

where the angular velocity of a field line f2p(\I/) and the rest-mass energy loading rate per 
unit magnetic flux k(^f) are also the integrals of motion. (These four integrals of motion 
for outflows are assumed to be positive in the following.) We call the toroidal field, and 
the poloidal field strength is denoted by B p . If the flux function ^(R, Z) is determined, the 
poloidal field strength is given by 

B 2 = (B R ) 2 + (B Z ) 2 , B R = --^, B z = -^, (4) 
p v ' v ; ' RdZ ' RdR ' v ; 

and the toroidal field can be expressed in terms of the poloidal quantities. Further, the 
Lorentz factor 7 equal to E k corresponds to bulk plasma motion involving both poloidal and 
rotational components. The angular velocity of plasma is denoted by Q in L k . The poloidal 
4-velocity u p = 71^ is used instead of the 3-velocity v p . Then, the integral of motion k is 
given by k = pu p /B p , where p is the proper mass density of plasma. 

Now we can define the relativistic Alfven Mach number M by the equation 

M 2 = Anpu 2 p /B 2 = 4irku p /B p , (5) 



which is used to give the toroidal field as 



-7- 



and x = = R/Rl. Because 7 and RQ can be also represented by the four integrals of 
motion and M 2 , the condition 7 2 (1 — R 2 Vt 2 ) — u 2 = 1 for normalization of 4-velocity leads 
to the poloidal wind equation 



1 + u 2 p )(l -x 2 - M 2 ) 2 = e 2 (l-x 2 - 2M 2 ) + (e 2 - M 4 , 



(7) 



where e = E — LQ F = 7(1 — R 2 QQ F ) is denoted by the plasma parts of the conserved energy 
and angular momentum. The angular velocity of plasma Q is given by 

_ (n F L/E)M 2 -x 2 (l-n F L/E) 

x(M 2 -l + n F L/E) ■ [) 

Let us rewrite equation (5) to give the poloidal 4-velocity u p in equation (7) as follows, 

u p = B p M 2 /4nk . (9) 

Then, by assuming a specific flux function Z) (giving the field strength B p ), one may 
solve the algebraic equation (7) for M 2 to determine the evolution along a fixed field line. In 
general, however, such solutions giving M 2 as a function of x and ^ would not be trans-fast- 
magnetosonic: If u p becomes equal to the relativistic fast-magnetosonic wave speed at a point 
on a field line of \1/ = constant, the partial derivative dM 2 /dx diverges there. (Otherwise, 
the flow remains sub-fast-magnetosonic or super-fast-magnetosonic everywhere on a field 
line.) It becomes necessary to find a special class of solutions satisfying the critical condition 
to keep a finite acceleration of u p at the fast-magnetosonic point (see, e.g., Takahashi & 
Shibata 1998). In this paper we do not adhere to such an analysis of the critical condition 
relating the four integrals of motion with B p at the fast-magnetosonic point. We rather turn 
our attention to a different useful information concerning plasma acceleration of trans-fast- 
magnetosonic flows, which is derived from the poloidal wind equation without assuming any 
details of a flux function ^/(i?, Z). 

For this purpose we introduce the poloidal electric to toroidal magnetic field amplitudes 
ratio £ as follows, 

^Jl = ^' (10) 

where we have used the relation E p = xB p between the poloidal field amplitudes. If equation 
(10) by the help of equations (6) and (9) is used, the poloidal velocity u p (and the Lorentz 
factor 7) is given by M 2 and £ (instead of B p ). Then, the poloidal wind equation (7) can be 
reduced to the quadratic equation for M 2 



AM' - 2BM 2 + C = , 



(11) 
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where the coefficients A, B and C are given by 



A = E>-l-e(E-^ -i^E, (12) 

B = x 2 + e 2 - 1 , (13) 
C = -(x 2 - 1)5 . (14) 

If the ratio £ is given as a smooth function of x under a constant the evolution of M 2 
along a field line can be derived from equation (11). Then, we can determine the toroidal 
field B^ from equation (6) and also the poloidal field B p from the relation B p = ^\B,p\/x. 
As will be briefly explained later, this field amplitude B p obtained through this procedure 
automatically satisfies the critical condition at the fast-magnetosonic point. This means 
that the quadratic equation (11) is not singular even at the fast-magnetosonic point. Hence, 
our necessary task in this section is just to solve explicitly equation (11) and to discuss the 
parametric range of £ for allowing the existence of highly relativistic trans-fast-magnetosonic 
solutions. 

After some algebraic manipulation, it is easy to show that equation (11) gives the 
solution of the form 

r 2 - 1 

M 2 = L —== , (15) 

[E-(m F /x")]^f/B-l 

where the function / depending on £ and x is defined by / = (1 — £, 2 )x 2 + £ 2 . It is clear that 
the function B involved in equation (15) cannot be negative, because e is estimated to be 

(i - ggtf 

- l-RHP ' [ ' 

Then, we must require the condition / > 0, under which we obtain the allowed range of £ 
for outflows propagating at a radius x > 1 beyond the light cylinder as follows, 

e<^- v in) 

A more stringent constraint is obtained if we consider M 2 at the asymptotic radius x — > oo, 
where we have 

M 2 1 

—2 -> / • 18 

x 2 

This leads to the result that if the inequality £ 2 > £ 2 = 1 — (1/-E" 2 ) always holds on a field 
line, M 2 becomes infinitely large at the radius x = x c given by 



Xr 



(19) 
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To avoid the divergence of M 2 at a finite radius, the field lines with £ > £ c should become a 
closed loop or be asymptotically cylindrical within the radius x = x c . Note that outside the 
light surface the closed loop configuration should be forbidden because of the inertia effect of 
plasma. So the flow streaming along the closed loop solution should transit to another open 
field solution by making a MHD shock. Then, we can expect that the asymptotically cylin- 
drical configuration is a more plausible one. (For the closed loop aligned flows without the 
shock formation, the ideal MHD approximation may be broken at least near the equatorial 
plane, because non-ideal MHD flows can stream across the magnetic field lines.) 

The Alfven point on a field line is known to be present at M 2 = M| w = 1 — x 2 , from 
which under the condition E > LQp the solution (15) gives the position x 2 = x\ = LQp/E 
and the Alfven Mach number M 2 = M\ = M% N {xx) = e/E. In general, the differential 
form of the poloidal equation contains a singular term (see, e.g., Takahashi & Shibata 1998). 
However, with the above condition, apparently the sub-Alfvenic outflows given by the wind 
solution (15) can smoothly pass through the Alfven point and propagate to the super- Alfvenic 
region x > x\. After passing through the light cylinder x — xl = 1, which is not a singular 
point for physical flows, the super- Alfvenic outflows would reach distant regions. Of course, 
M 2 remains finite also at the light cylinder, where we have 

Ml s M^ L ) = J^L, w = l±«!±i)f! , (20) 
L v ; 1 - w 2eE K J 

The Alfven Mach number at the light cylinder becomes larger if compared with the value 
Ma at the Alfven point. 

Now let us give the field amplitude B p , using M 2 obtained by equation (15). From the 
calculation of and the definition of £ we find the simple relation 

Ank£ [b . x 

B p = ■ 21 

p M 2 y / y ' 

This allows us to estimate the field amplitude B p at the fast-magnetosonic point x = xf 
with M 2 = Mp = Mp W (x F ), where we have 



= ^ (22) 



for the Alfven Mach number related to the fast magnetosonic wave speed. From equation 
(15) it is easy to find that the relation M 2 = M| is satisfied at the point determined by the 
equation 

E-^f) 2 f = X 'B. (23) 
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We must remark that the light cylinder radius x — 1 satisfying equation (23) is not the 
fast-magnetosonic point. In spite of the existence of such a redundant solution, we use this 
concise form (23) for mathematical simplicity. In fact, equation (23) is helpful for calculating 
the partial differentiation of equation (21) with respect to x along a fixed field line at the 
fast-magnetosonic point, and we can straightforwardly check that the terms involving d^/dx 
automatically cancel out in the derivative dB p /dx at this critical point, and the critical 
condition for the value of the partial derivative of B p holds, irrespective of any choice of 
£. (If the derivative <9£/cfe diverges at the fast-magnetosonic point, the critical condition 
for dBp/dx is not satisfied. This means that a smooth change of £, namely B^, at the 
critical point is the essential requirement for a trans-magnetosonic flow.) Thus, we can 
obtain a trans-fast magnetosonic flow solution without the regularity condition at the fast 
magnetosonic point. The behavior of the acceleration in the outgoing flow is determined by 
£ = £{ x )i an d the asymptotic feature depends on the value of £; that is, £ > £ c or £ < £ c . The 
typical solutions for outgoing trans-fast magnetosonic flows are demonstrated in Figures la 
and lb, for both cases of £ < £ c and £ > £ c . The outflows can start from the plasma source 
(x <C 1), and after passing through the Alfven point (A) and the fast magnetosonic point 
(F) they reaches distant regions. In the case of £ < £ c , the flow reaches x 3> 1 region with a 
finite Mach number. This is the standard picture discussed by lots of previous wind models. 
On the other hand, in the case of £ > £ c the flow confines within x < x c . We can expect 
that asymptotically the magnetic field line becomes cylindrical and the flow streams toward 
Z-direction. In the distant Z-region (ZQp ^> 1), the flow becomes to have a very high 
Alfven Mach number. 



2.2. Acceleration and Energy conversion 

Though we have constructed the basic formulae for discussing the evolution of M 2 , in 
this paper we are particularly interested in highly relativistic acceleration of bulk motion 
through a conversion of magnetic to kinetic energy. In a sub-Alfvenic region the outflows are 
expected to be injected under a magnetic-energy dominated state (E ~ LQp) with a very 
large value of E. This means that we can analyze the evolution of M 2 in more details, using 
the approximations 

r - i o, 

< 1 , E > 1 . (24) 



e E — m F 



E E 

For example, equation (23) claims that the value of x 2 ? at the fast-magnetosonic point be- 
comes very large and is approximately given by E 2 / 3 fp, in which the exact form of / = f{x) 
is necessary because £p = £(#f) may be very close to unity. Hence, we can write the position 
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of the fast-magnetosonic point as follows, 

£2£2/3 

= 1 + - 1)EW ■ (25) 

We note that this critical point exists only if £| > 1 — (1/E 2 ^ 3 ), while the allowed range is 
1 < xp < oo, depending crucially on the difference £p — 1. On the other hand, for the value 
of Mp at the fast-magnetosonic point we obtain 

f = ^«i- («•) 

To show the physical implication of equation (26), we consider the specific magnetic 
energy E m at a radius x ^> 1 beyond the light cylinder. From equations (2) and (6) we 
obtain approximately 

Em = T+W' (27) 

where we have introduced the quantity M 2 defined by M 2 = M 2 /x 2 . Then the ratio of the 
specific kinetic energy E k = E — E m to E m is given by 

|^ = M 2 , (28) 

which implies that M 2 is an indicator of the conversion of the huge magnetic energy into the 
kinetic energy for outflows propagating to a radius x 1. (The term M 2 corresponds to 
the inverse of the magnetization parameter a in some wind models; see, e.g, Michel (1969) 
and Camenzind (1986) for the radial wind, and Begelman & Li (1994) and Takahashi & 
Shibata (1998) for non-radial winds.) At the fast-magnetosonic point, however, we obtain 
from equation (26) 

^)r^ <<u (29) 

which is a generic result independent of a field configuration ^(R, Z) of highly relativistic 
outflows. It is clear that the energy conversion is still inefficient at the fast-magnetosonic 
point. 

Then, we consider plasma acceleration of trans-fast-magnetosonic outflows propagating 
at a radius x ^> 1 beyond the critical point, where we have approximately 

M 2 = 1 . (30) 

y/o- - e)E 2 + (eE 2 /x 2 ) - 1 
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As was previously mentioned, M 2 becomes infinitely large at x = x c if £ > £ c . For highly 
relativistic outflows this critical radius is approximately given by 

X * = 1 + (£2 _ l)E 2 ■ ( 31 ) 

A possible configuration of a field line given by the parameter £ larger than £ c would be 
asymptotically cylindrical, because it can be confined within a radius smaller than x c even 
in the limit Z — > oo. Unless £ is very close to unity, the value of x c is of order of unity, 
and such a cylindrical jet becomes very narrow. In the limit x — > x c , the energy ratio M 2 is 
estimated to be 

M ' = 1 + (£2 _ 1)E 2 X 1 _ {x / Xc) > ( 32 ) 

The narrow jet corresponding to £ c of order of unity can be kinetic-energy dominated only 
if the asymptotic radius x is very close to x C) i.e., 1 — (x/x c ) = 0(1/ E 2 ). This is the fine- 
tuning problem required for the cylindrical field to realize highly relativistic acceleration of 
bulk motion. On the other hand for a cylindrical field line with the value of £ such that 
|£ 2 — 1| = 0(1/E 2 ) the maximum radius x may become of order of E 2 , and from equation 
(30) we find that at least a rough equipartition ~ E m will be realized at the scale of 
x ~ E without the fine-tuning of x — > x c . In particular, if £ < £ c , namely, £ 2 < 1 — (l/E 2 ), 
the field line may be extending to an infinite radius x — > oo, where M 2 becomes equal to 
(Ey/l — £ 2 — l) -1 . Such a field line configuration may be asymptotically paraboloidal or 
conical, and the efficient energy conversion into the state Ej, > E rn becomes possible at 
x > E for field lines with 1 — £ 2 = 0(1/E 2 ). This is another fine-tuning problem required 
for the poloidal electric and toroidal magnetic field amplitudes. The precise fine-tuning of 
£ 2 = 1 — (l/E 2 ) means that the magnetic energy E m can be completely transported into the 
kinetic energy E k as outflows propagate to an infinite radius. 

To claim that the kinetic energy can asymptotically become larger than the magnetic 
energy for injection of magnetic-energy dominated outflows with very large E, we must solve 
the fine-tuning problem such that 1 — (x c /x) = 0(1/E 2 ) or |£ 2 — 1| = 0(1/E 2 ) as a result 
of MHD interaction described by the Grad-Shafranov equation. In this paper we consider a 
jet ejection with a very small opening angle such that R/Z < l/E, and we study the spatial 
variation of £ to show the dynamical fine-tuning of £ in jet flows. 



3. THE APPROXIMATED GRAD-SHAFRANOV EQUATION 



The asymptotic analysis of relativistic outflows has been developed in previous works 
(Li, Chiueh & Begelman 1992; Appl & Camenzind 1993a,b; Eichler 1993; Begelman & Li 
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1994; Tomimatsu 1994), and the logarithmic dependence of the asymptotic structure on x 
has been pointed out. However, the transition from magnetic-energy dominated state into 
kinetic-energy dominated one is an important unsolved problem which is beyond the usual 
scheme of the asymptotic analysis based on the naive approximation x ^> 1 in the poloidal 
wind and Grad-Shafranov equations. In our approach presented in the previous section, such 
an approximation corresponds to 

/~(1-£V»1, (33) 

for which we have 

M 2 = -— (34) 

from equation (30). (This should be the case of £ < £ c .) The key point missed in this 
calculation is that the fine-tuning of 1 — £ 2 = 0(1/ E 2 ) to realize the energy equipartition 
(M 2 ~ 1) may occur at a radius x in the range 1 < i < E. Then, to discuss the energy 
conversion in the range 1 <C x < E, which is called the "intermediate" range of x in this 
paper, we must analyze the evolution of M 2 without assuming (1 — ^ 2 )x 2 to be very large. 
The existence of the intermediate range of a; is a main feature of highly relativistic outflows 
with a very large specific energy (E 3> 1). Recalling that we obtain M 2 ~ e/E at the light 
cylinder surface x — 1, we expect M 2 to increase from a sub-fast-magnetosonic value in the 
range 1/E <C M 2 <C 1/ E 2 ^ to a rough equipartition value of M 2 ~ 1 as outflows propagate 
in the intermediate region. The fast-magnetosonic point at which we have the value of 
M 2 = l/E 2 / 3 can be involved in this region. The usual asymptotic analysis becomes valid 
only in the region x ^> E, where M 2 may increase logarithmically with x, and our purpose 
here is to give the more precise treatment valid both in the intermediate and asymptotic 
regions. 

The Grad-Shafranov equation can be written in the form 

M 2 + x 2 -l /M 2 + x 2 -l \ 

-V W = Sl + s 2 , (35) 



8vr 2 V R 2 

if no gravity is included (see, e.g., Tomimatsu 1994). The source terms s\ and s 2 are given 
and 

s2 = .Ml±^l^ _ *^ (M2 + ^ - 1 + ■ (37) 

Here we regard M 2 , d^/dR = ^ R and d^/dZ = ^ z in equation (35) as functions of the 
variables R and \I> ', instead of R and Z. Further, by virtue of the introduction of £, from 
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equations (15) and (21) we find the relation 



M 2 + x 2 - 1 = Ank IE 



LSI, 



x c 



B, p 



(38) 



which is substituted into the left-hand side of equation (35). Then, after some manipulation 
we obtain 



J__d_ 



ek 2 

l + q 2 



E — 



F 



X, 



Si + S 2 , 



(39) 



where q = — ^ z/^r means the slope of a poloidal magnetic field line, and the modified 
source term s\ is 



si = - 



d_ 

mi 



r 



2 f 2 



(40) 



Now we derive the approximated form of the Grad-Shafranov equation valid for highly 
relativistic outflows (with E ^> e > 1) propagating in the intermediate and asymptotic 
regions, where we obtain 

2 



l>l-^+^r~— l + ^r- . 

x 2 E 2 \ M 2 



(41) 



as a result of the approximation M 2 + x 2 — 1 ~ [M 2 + l)x 2 in equations (38) and (21). It is 
interesting to note that the source terms are reduced to the compact form 



2 , , ,2 9 ( k 2 



§! + s 2 — (1 + M )VL F — 



\ Q 2 F M 4 

To see clearly the fine-tuning of £, we rewrite it as follows 

,2 i P 
c ~ 1 



where we define 







M 2 ) 



x 2 



(42) 



(43) 



(44) 



The renormalized variable x = x/E is also useful to discuss the evolution of M 2 in the 
regions considered here. 

The observations of AGN jets (and possibly GRB jets) show the existence of highly 
collimated outflows with very narrow opening angles, for which we can assume that q 2 <^ 1 
in equation (39). Hereafter, we consider only the case of q < l/E to give 

AnkE 



m R ~ r 2 b p ~ 



n 2 F (i + M 2 ) 



(45) 
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This allows us to eliminate from the left-hand side of equation (39), in which we also 
have [E — (LVL F /x)} 2 ~ E 2 — (2/x 2 ). Then, we arrive at the final form of the reduced 
Grad-Shafranov equation 

df3_ _^E_d_ ( k 2 \ A_ 

OR = k \tf F M* J + ~x 2 ' ( j 

If equation (44) is substituted into equation (46), we can obtain the partial differential 
equation of first order for M 2 as follows, 

R(M 2 + 1) DM 2 b dM 2 _ a 1 

Af6 dR + M6 ~ M4 x 2 ' ^ ' 

where 

A*E3 (k \ 



and 

In the intermediate region (x < 1) the role of the term 1/x 2 in equation (47) becomes quite 
important, while it has been completely neglected in the asymptotic analysis corresponding to 
the case of x 3> 1. If a solution M 2 (R : \l>) is derived from equation (47), it is straightforward 
to find the flux function \l/(-R, Z) from equation (45). In the next section we will present the 
interesting example of a jet solution to discuss the energy conversion and the change of field 
configuration in the intermediate and asymptotic regions. 



4. FIELD STRUCTURE OF JET FLOWS 

Equation (47) is valid only for jets with small opening angles propagating beyond the 
light cylinder surface. Considering the limit * — > for such outflows, we expect the jet 
ejection to occur near the polar axis RQ F — > 0, where the toroidal filed should decrease 
in proportion to R. Noting that the small flux function \1/ is also proportional to R 2 , we 
estimate to be E ~ E m ~ ty/k in the injection region near the polar axis. Then, to obtain 
a large specific energy E per one particle, the rest-mass energy loading rate k per unit 
magnetic flux should also become small in proportion to This boundary condition near 
the polar axis motivates us to consider the case such that E and £l F are independent of ^, 
while k = k (ty/^o), where k and ^ are constants and should be given by the boundary 
conditions at some foot point in a plasma source. Then, we obtain 
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where a is a dimensionless constant. 

Further, let us recall that M 2 ~ 1/E near the light cylinder surface, Then, we expect 
M 2 to be very small in proportion to x in the region l/£ C x C 1, from which equation 
(45) leads to R^r ~ a$> . This means that for a = 2 we have \I/ oc R 2 , assuring the smooth 
matching to the inner solution valid near the polar region. (Of course, we cannot give the 
condition a = 2 to be necessary, because our analysis is limited to the outer solution valid 
in the range RQp ^> 1.) Fortunately, for the model with a = 2 and b = 2\1/, we can give the 
general solution of an analytical form for equation (47) as follows, 

2f 2 (l+M 2 ) , /25 2 \ n 

In U- + 1 + D l , (51) 



2x 2 + M 2 \M 2 
where D\ is an arbitrary function of D 2 defined by 

D 2 = — ( J- + -L} . (52) 
\M 2 2x 2 J 1 ; 

Applying equation (45) to this general solution, we find that D 1 and D 2 are arbitrary func- 
tions of Z, which should be determined by the additional boundary conditions for jet flows. 
We can study an essential feature of plasma acceleration in jet propagation from this simple 
model. 

First let us discuss the results obtained from equation (51). Considering that M 2 ~ 1/E 
near the light cylinder x L = 1/E, we claim that x 2 /M 2 ~ x <^ 1 in the range 1/E < x <^ 1, 
from which equation (51) leads to 

2x 2 (l — x 2 /M A ) ~ Di(Z) . (53) 

for all Z. This equation means that the absolute value of D\(Z) should be chosen to be 
at most of order of 1/E 2 for all Z. In fact, if \Di\ ^> 1/E 2 for some Z, we can easily see 
that equation (53) breaks down in the range l/£ 2 C i 2 < l-^il- Hence, we can discuss 
the increase of M 2 ~ Ek/E m in the range x ^ 1/E under the choice of D\ = 0, namely, 
according to the following equation 

2^ + l) ^ A 



2P + M 2 \M 2 

It is now clear that we have M 2 ~ x in the range l/£ « x « 1, and the outflows can 
pass through the fast-magnetosonic point (i.e., M 2 = 1/E 2 / 3 ) located on the radius RUf — 
E 1 / 3 . If the outflows can arrive at the radius x ~ 1.4, the equipartition E k = E m between 
magnetic and kinetic energies is realized. In the asymptotic region a; > 1 we can confirm 
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the logarithmic increase of M 2 given by M 2 ~ ln(2x 2 /M 2 ) toward the complete conversion 
of magnetic to kinetic energy. The numerical solution of equation (54) is shown in Figure 
2a. The corresponding poloidal velocity u p (x) and Lorentz factor 7(x) are also shown in 
Figure 2b, where the Lorentz factor is given by 7 = E — u p /(C,M 2 ). In this case, we obtain 
< 1 anywhere for the ratio of the poloidal electric to toroidal electric magnetic field 
amplitude (see Fig. 2c). Note that the Lorentz factor includes both the poloidal and toroidal 
motion. In Fig. 2b, the difference between -u p -value and 7- value near the light cylinder is due 
to the dominated toroidal motion of the plasma. (Just around the light cylinder x ~ x L , 
the properties of the flow shown in Fig. 2 may be incorrect because of our approximations, 
but we can expect correct features at least near and outside the fast magnetosonic point, 

x > Xp. ) 

Next, we discuss the field configuration given by equation (52). The jet flows may be 
confined by an external pressure (see, e.g., Li 1993; Begelman & Li 1994; Fendt 1997). If 
the shape ZQp = H(x) of the last flux surface ^ = ^0 is determined by the outer boundary 
condition, the function DiiZi) in equation (52) is fixed. For example, let us consider the 
radial last flux surface at an angle R/Z = #0 with the pole axis direction. (Because the 
solution can be applied only to jets with small opening angles, we must require that $o is 
at most of order of 1/E.) Using the function M 2 = M 2 (x) derived from equation (54), we 
obtain ^ 

^^iw^ (55) 

along a flux function fy(R, Z) = constant, where z = ZQp/E. Note that for the value of 
M 2 (6 Q z) in equation (55) we use the function M 2 = M 2 (x), where the variable x should be 
replaced to 9 z. Then, in the range 9 ZQ F /E <C 1 of Z, equations (52) and (55) leads to 
the conical field configuration 

(56) 

This shape of field lines changes as Z increases, and in the asymptotic region we find the 
leading behavior such that 

6f)Z^lp ( R^lp^ 



E ~ V E • ■ < 57 > 

which represents a paraboloidal collimation of the inner field lines for \1/ < \l/o- Figure 3 shows 
the shape of the magnetic field lines, which can be solved numerically by using equations 
(55) and (52) with the function M 2 = M 2 (x) ( or M 2 = M 2 (6 z) ) derived from equation 
(54). 

From the solution (54) we note that the ratio M 2 of kinetic to magnetic energy is 
completely determined by the propagation radius x = Rflp/E, irrespective of the field 
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configuration given by equation (52). If the last flux surface is asymptotically cylindrical 
with the outermost radius x , giving the function D 2 (Z) approaching a constant value in the 
limit Z — > oo, the maximum value of M 2 in the cylindrical jets is crucially dependent on x . 
For narrower jets with x$ < 1 the magnetic energy is still dominant even in the asymptotic 
region. 

The outer flux surface corresponds to a larger radius x if compared with a fixed vertical 
distance Z. Thus, we can claim also that the energy conversion along outer flux surfaces of 
larger \1/ becomes more efficient. This tendency of the efficient energy conversion is due to 
the condition for the rest-mass energy loading rate per unit magnetic flux such that k oc 
In fact, if we consider the case that E, Qf and k are independent of we find a different 
dependence of M 2 on x from the general solution for equation (47) with a = and a constant 
b, which is given by 

i V i 



and 



M 2 J x 



~2 r2 I M 2 F 1 - M 2 - 1 



l/Fi 



exp — = x 2 F 2 5 . (59) 

V b ) \M 2 Fi + M 2 + l/ 

From equation (45) we can check that F x and F 2 are arbitrary function of Z . If we assume the 
behavior of these functions in the limit Z — > oo such that Fi ~ 1 + (1/lnZ) and F 2 remains 
finite, we can confirm the increase of M 2 in logarithmic scale of Z and the paraboloidal 
shape of field lines similar with equation (57) in the asymptotic region. However, it is clear 
from equation (58) that M 2 decreases as x increases under a fixed Z. The energy conversion 
becomes less efficient on outer flux surfaces. Further we note that the solution (58) fails to 
give a real value of M 2 in the range i<1 under a fixed Z. Such a region should be covered 
by inner flux surfaces corresponding to small and the previous solution (54) with k oc ^ 
should be used there. One may consider a more realistic dependence of k on Then, a 
smooth change of M 2 from equation (54) to equation (58) will be allowed as ^ increases in 
the range < * < 

It is sure that according to a choice of the integrals of motion as functions of * we can 
obtain various evolutionary models of energy conversion in jet flows, which may show more 
complicated behaviors different from the simplest solution (54). However, we would like to 
emphasize the key result obtained here that a rough equipartition between magnetic and 
kinetic energies is realized at the radius R ~ RlE far beyond the light cylinder radius i?L, 
and the subsequent logarithmic increase of kinetic energy goes on at larger radii R ^> R-^E. 
This potentiality of MHD acceleration will be robust at least under the boundary condition 
for inner flux surfaces ^ — > such that the rest-mass energy loading rate is given by k oc ^ 
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to keep the total specific energy E very large. If so, we can claim that there exists the 
critical scale given by R = RlE for conversion of Poynting flux injected in jet flows, which 
is important for discussing a prompt emission of radiation owing to dissipation of kinetic 
energy of bulk motion. 

5. SUMMARY AND DISCUSSION 

We have presented a parametric representation of Alfven Mach number M by the ratio 
£ of poloidal electric to toroidal magnetic field strength, to avoid the troublesome analysis of 
the critical condition at the fast-magnetosonic point. Then, from the poloidal wind equation 
we have easily derived trans-fast-magnetosonic solutions including the parameter £ defined 
as a smooth function along a field line. If the parametric function £ is asymptotically larger 
than the critical value £ c , the flux surface have been shown to be confined within a finite 
radius x = x c . Otherwise, the dynamical fine-tuning of £ have been required for acceleration 
to highly relativistic bulk speeds. To determine the parametric function £, we have also 
given the approximated form of the Grad-Shafranov equation applied to jets ejected with a 
very large total specific energy E and confined within a very small opening angle of order of 
l/E. 

Acceleration of outflows in generic models with the integrals of motion E, Q F and k 
variously dependent on \P has not been analyzed in this paper. The comparison of the 
two models given by equations (54) and (58) suggests a variability of acceleration efficiency 
depending on various choices of the integrals of motion. Nevertheless, the former solution 
(54) studied in the previous section will be interesting as a typical model revealing the 
high potentiality of MHD acceleration, which was also discussed by Okamoto (2002) in 
relation to pulsar winds (see also Michel (1969); Begelman & Li (1994)). By virtue of the 
closed-form expression (54) of the model, we can clearly understand the following evolution 
of jet flows in the intermediate and asymptotic regions far beyond the light cylinder, if 
the jet radius R extends to an infinite distance: (1) The magnetic-energy dominated and 
sub-fast-magnetosonic outflows pass through the light cylinder surface (R — Rl) with the 
Alfven Mach number such that M 2 ~ l/E. (2) Then, the energy ratio Ek/E m increases 
to l/E 2 / 3 at the fast-magnetosonic point corresponding to the radius R ~ RlE 1 / 3 . In 
the intermediate region (i?L <C R < RlE) the outflows can smoothly become super-fast- 
magnetosonic. (3) Importantly, the model claims the realization of rough equipartition 
between kinetic and magnetic energies at the radius of order of RlE. (4) The further 
energy conversion toward a kinetic-energy dominated state in logarithmic scales of R in the 
asymptotic region (R 3> RlE) is also confirmed. Figure 4 is the summary of this jet solution. 
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The full conversion of magnetic energy into kinetic one means that the asymptotic 
Lorentz factor of bulk motion becomes equal to the total specific energy E of outflows 
injected near the central source. Hence, for kinetic-energy dominated jets observed with a 
huge bulk Lorentz factor 7, we can expect the rough equipartition ~ E m to occur at the 
jet radius Rj Ct = -Rl7 [ Note that the value of 7 (~ Ek) is roughly same order with E ]. 
Though we have considered ideal MHD flows in this paper, the observed jet activity, such as 
a prompt emission of radiation and a ultra relativistic acceleration of electrons, should be due 
to dissipation of the power of bulk motion, for example, through formation of shocks. Then, 
the interesting high-energy phenomena of jets will be observed, only after the kinetic energy 
of bulk motion begins to dominate, namely, the jet radius extends to this critical radius Rj et . 
The shock formed in the energy equipartition region would be distinct from shocks formed in 
the kinetically dominated asymptotic region by observations. For the kinetically dominated 
flows, the compression ratio behind the shock is much higher, and flatter synchrotron spectra, 
higher emissivities, etc. would be observed (Begelman & Li 1994). 

In particular, for AGN jets with 7 ~ 10 (see, e.g., Ghisellini et al. 1993), we can estimate 
the critical radius to be i?j e t ~ 10/?l- When magnetic fluxes for the jet connect to a rotating 
geometrically thin disk around a black hole, we can regard Qp as Q^(r), where f2i<:( r ) is the 
angular velocity for circular equatorial orbit in the Kerr metric, which corresponds to the 
Keplerian angular velocity in the Newtonian case. If the foot points of most magnetic fluxes 
distribute near the inner part of the disk r ~ r ms , the angular velocities of the magnetic field 
lines are roughly ~ ^K<7ms) (see, e.g., Camenzind & Krockenberger 1992). Then, the 

light cylinder radius R L = l/Q F for field lines threading a rotating black hole and the disk's 
inner part is possibly equal to several Schwarzschild radii. Hence, our conjecture from the 
jet model obtained here is that the active region in AGN jets, far from the central source, 
has the radius as large as a hundred of Schwarzschild radii, corresponding to 0.01 pc for a 
10 9 M Q black hole. This result is consistent with observational results of blazars (Kataoka 
et al. 2001) and of the M87 jet (Junor, Biretta, & Livio 1999). (For GRB jets expected to 
have higher Lorentz factors 7 > 100 of bulk motion, the size of the emission region will be 
larger than a thousand of Schwarzschild radii.) The dependence of the radius of the emission 
region on the bulk Lorentz factor is an interesting problem to be checked by observations of 
AGN jets. 

Finally let us discuss the vertical distance Z- )et from the central source, which gives the 
critical radius R- iCt on a field line. We have postulated jet flows confined within a small 
opening angle of order of 1/E and derived the solution valid only in the region Z > ER 
from the Grad-Shafranov equation. We can obtain the field configuration according to 
equation (52) under a suitable boundary condition on the outermost flux surface. As an 
example, in the previous section, we have considered the case with the radial outermost 
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flux surface given by Z = R/9 . If the opening angle 9q is equal to I/7, we can claim that 
the kinetic energy of bulk motion begins to dominate at the vertical distance Z Jct larger 
than 7-Rjet = 1 2 Rl, far from the central source. (The value of Z Jct becomes larger on inner 
flux surfaces corresponding to smaller ^.) This will predict a sub-parsec distance to the 
active region in an AGN jet from a 10 9 M Q black hole. We expect that at this distant 
region (Z > Zj Ct ) internal shocks form and make the active region in the jet. The internal 
shock scenario is generally thought for gamma-ray bursts, radio-loud quasars and blazars 
(see, e.g., Mezaros & Rees 1993; Spada et al. 2001; Ghisellini & Celotti 2002). In this 
scenario, it is assumed that the active region powered by collisions among different part of 
jet itself, moving at different bulk Lorentz factors. Although we consider stationary flows 
in the magnetosphere, we can also expect some kinds of plasma instability (e.g., the screw 
instability of plasma in black hole magnetosphere discussed by Tomimatsu, Matsuoka & 
Takahashi (2001)) to make blobs of plasma in the jet. If the cylindrical collimation of the 
outermost flux surface develops in the range R < i?j Ct , however, the vertical scale Zj et should 
be much larger. Because we have not discussed physical mechanisms (possibly due to external 
matter) of the confinement of the magnetic flux, the boundary condition on the outermost 
flux surface to determine the field configuration remains arbitrary. This is a problem beyond 
the scope of the present paper, which is important for giving a more definite estimation of 
the vertical distance from the central source to the active region in AGN jets. 

This research was partly supported by Scientific Research Grant (C2) 14540253 of Japan 
Society for Promotion of Sciences. 
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Fig. 1. — Examples of outgoing flow solutions (thick curves) given by equation (15) with 
(a) f = 0.8 and (b) f = 1.1, where E = 10.0 and x A = 0.8. The curves of the Mach 
numbers corresponding to the Alfven wave speed and the fast-magnetosonic wave speed, 
namely, M 2 = M| w and M 2 = M| w , are also plotted, where M| w = 1 - x 2 and M| w = 
1 — x 2 + (x 2 /C, 2 ). The crossings of these curves with the flow solution labeled by "A" and 
"F" are the Alfven and fast-magnetosonic points, respectively. 
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logx 

Fig. 2. — A solution of the approximated Grad-Shafranov equation given by equation (47) 
with E = 10.0. (a) The Alfven Mach number with two limiting curves (M 2 ~ x for x L <C 
x C 1 and M 2 ~ ln(25 2 /M 2 ) for x 3> 1), (b) the poloidal velocity u p (x) (solid) and the 
Lorentz factor 7(3:) (dashed), and (c) the ratio of the poloidal electric to toroidal electric 
magnetic field amplitude This solution is valid for x > x L = 1/E. 
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Fig. 3. — The magnetic field configuration corresponding to the solution given by Fig. 2. It is 
assumed that the jet flow is confined within a cone of small opening angle 9 by an external 
pressure. For the inner field lined of \I> < a paraboloidal collimation is obtained. 
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Fig. 4. — Schematic picture of the jet solution given by Fig.3. The term 6q is the opening 
angle of the last magnetic flux surface. The initial magnetically dominated outflow E m 3> E k 
accelerates along an almost radial field line at the intermediate region x L < x < 1 (region 
I), and then a rough equitation E k ~ E m is realized around the radius x ~ 1 far beyond the 
light cylinder (region II). Subsequent logarithmic increase of E k goes on at large radii x ^> 1 
(region III) along a collimating field line. 



